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Given a list of complex numbers Ai, . . . , A n , when can it be the spectrum of a quantum channel, i.e., a 

completely positive trace preserving map? We provide an explicit solution for the n — 4 case and show 

that in general the characterization of the non-zero part of the spectrum can essentially be given in terms of 

its classical counterpart — the non-zero spectrum of a stochastic matrix. A detailed comparison between the 

__ classical and quantum case is given. We discuss applications of our findings in the analysis of time-series 

and correlation functions and provide a general characterization of the peripheral spectrum, i.e., the set of 

eigenvalues of modulus one. We show that while the peripheral eigen-system has the same structure for all 

k^ Schwarz maps, the constraints imposed on the rest of the spectrum change immediately if one departs from 

pL** complete positivity. 
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I. INTRODUCTION AND OUTLINE 

Quantum channels provide general input-output relations for quantum mechanical evolutions. 
If the input space equals the output space, we can assign a spectrum to each quantum channel. 
The spectral properties play a role wherever successive powers of the channel are applied. These 
may describe discrete time-evolution or they may appear in particular representations of quantum 
spin chains. Apart from the most basic facts which coincide with their 'classical' counterparts for 
stochastic matrices, not much is known about the spectral properties of completely positive maps. 



The present work aims at a deeper analysis, beyond Perron-Frobenius theory, and provides nec- 
essary, and in some cases also sufficient, conditions for a list of numbers to be the spectrum of 
a quantum channel. That is, we address the 'inverse eigenvalue problem' for completely positive 
maps — occasionally being content with plain (not complete) positivity. 

After having collected some basic results in Seclll] SecllHlprovides a rather exhaustive analysis of 
qubit-maps — a complete characterization of the spectra of (completely) positive qubit maps is given. 
In Sec|IV|the non-zero part of the spectrum of general (completely) positive maps is investigated and 
close analogies to the classical case of stochastic matrices are found. SecfV]provides a detailed study 
of the peripheral spectrum (i.e., eigenvalues whose modulus equals the spectral radius) and of the 
corresponding eigenvectors. Before outlining some applications in the analysis of time-series and 
correlation functions in Sec|VH|we compare the quantum and the classical world, i.e., the spectral 



properties of quantum channels and those of stochastic matrices in Sec VI 



II. BASIC SPECTRAL PROPERTIES OF QUANTUM CHANNELS 

We begin with fixing some notation and recalling basic results. Let T : A4d(C) — > .M<z(C) be 
a linear map on the space of complex valued d X d matrices. T is said to be positive if it maps the 
cone of positive semidefinite matrices into itself. If T is positive, then it is in particular Hermiticity 
preserving, i.e., T(Ay — T(A^) for all A e A^d(C). The adjoint map T* is defined by imposing 
tr [AT*(B)\ = tr [T(A)B] for all A,B £ M d (C). T is trace-preserving iff T* is unital, meaning 
that T* (1) = 1. A linear map on .Md(C) is completely positive iff T®idd is positive. A completely 
positive and trace-preserving linear map is referred to as quantum channel. 

As a vector space _Md(C) is isomorphic to C d so that T can be represented as a d 2 x d 2 matrix 
for which we will write T. Fixing a Hilbert-Schmidt orthonormal basis in .Md(C), i.e., a collection 



If T is 



of d 2 operators {E, £ M d (C)} with tr E\E 3 = 5 tj , we get T l3 = tr E\T(E. 

Hermiticity preserving, then T* is represented by the adjoint matrix T\ 

When analyzing the spectrum of a map we distinguish between the spectrum as a set, where each 
element only appears once, and as a multiset (unordered n-tuple), where elements appear according 
to their algebraic multiplicity. More precisely, the spectral set of T is the set of complex numbers A 
for which T— Aid is not invertible. The spectrum (as a multiset) is the list of roots of the characteristic 
polynomial det[T — Al] where each root appears with a multiplicity according to its degeneracy. 
We will write spec(T) if we refer to the multiset and {spec(T)} if the spectral set is meant. 

From Hermitian conjugation of the eigenvalue equation T(X) = XX we see that the eigenvalues 
of Hermiticity preserving maps are either real or come in complex conjugate pairs. This holds 
even if we take multiplicities into account. Since spec(T) = spec(T*), every trace-preserving 
map has an eigenvalue one. If the map is positive in addition, then its spectral radius is one, i.e., 
max{|A| |AGspec(T)} = l. 

A positive map is irreducible [EHK78 1 iff the spectral radius is a non-degenerate eigenvalue and 
the corresponding eigen-vector is positive definite. If in addition there is no other eigenvalue whose 
modulus matches the spectral radius then the map is called primitive [SPGWC09|. 

The following shows that a characterization of the spectra of (completely) positive maps is essen- 
tially independent of whether or not the maps are required to be trace-preserving: 

Proposition 1 (Similarity with trace-preserving maps) Let T : Md(C) -¥ Md(C) be a (com- 
pletely) positive map with spectral radius g > 0. There is a trace-preserving (completely) positive 
map T' on A^rf(C) such that spec(T) = q spec(T'). 

Proof. Define a map r e :l4 T(X) + etr [X] 1 and denote its spectral radius by Q e , Following 
BEHK78I T e is irreducible for all e > and for e — > we have that spec(T c ) — > spec(T) continu- 
ously. Irreducibility implies the existence of a positive definite P = T <l (P)q~ 1 > 0. Hence the map 



T' € {X) := e -ip-i/2i;(pi/2l-pi/2)p-i/2 satisfies spec(T e ) = spec(T^ £ and is (completely) 
positive and unital. Since these maps form a compact set, the limit T' = lim e _j.o (T'* ) exists and has 
all desired properties. |— I 

In the following we will consider trace-preserving maps only. 

III. THE n = 4 CASE - QUBIT CHANNELS 

We continue the discussion with completely positive qubit maps T : A^C) -» A^C) which 
allow for an exhaustive analysis. Some of the results derived in this section will be the basis for 
the general characterization of the non-zero part of the spectrum in Sec jIV A| We will first recall 
and derive some preparatory results and then provide an explicit solution of the inverse eigenvalue 
problem in Thmfl] 

Choosing normalized Pauli matrices as operator basis (with cto = 1) we can represent T as a 
4x4 matrix T^ := tr [<jjT(cr,-)] /2. If T is Hermiticity preserving, then T is real, and if T is 
trace preserving, then (Ti j) — (1,0, 0, 0). Since the maps of our interest satisfy both, they can be 
represented as 

where A is a real 3x3 matrix and ti € I 3 . The spectrum of T is thus the union 

spec(T) = {l}Uspec(A). (2) 

The Jamiolkowski state corresponding to T is given by r = j ^ . . Tij<Ji (g> aj so that v describes 
its reduced density operator and A its correlations. The (generally cumbersome) conditions for 
T to correspond to a completely positive map have to be read off from t > (cf. |RSW02|). 
The parametrization allows, however, for a nice geometric interpretation: parameterizing a density 
operator via p — (1 + J2k=i x kOk)l% i-e., in terms of a vector x £ M 3 within the Bloch ball 
| |x| I2 < 1, the action of T as parametrized in Eq.(fT]i is a simple affine transformation 

xt-*v + Ax. (3) 

From here conditions for T to be positive are readily derived as the vector has to stay within the 
Bloch ball. As the center of the Bloch ball corresponds to the maximally mixed state, v = holds iff 
T is unital. In this case positivity is equivalent to 1 1 A| |oo < 1, The following is a useful proposition 
which allows us to restrict ourselves to the unital case: 

Proposition 2 (Reduction to unital maps) Let T : A^2(C) — >• A^C) be a trace-preserving lin- 
ear map with Ay := tr [crjT(<Tj)] /2 the lower-right submatrix (i.e., i,j = 1, 2,3) of the matrix 
representation T. Then the map defined byT' :— lffiAis unital, trace-preserving and has the same 
spectrum as T. Moreover, it is (completely) positive ifT is. 

Proof. By construction T" is trace-preserving and unital and has the same spectrum as T, so it 
remains to show that (complete) positivity is preserved when going from T to T". To this end 
we use that D := diag(l, — 1, — 1, —1) is the matrix representation of time-reversal, i.e., matrix 
transposition in some basis. A look at the Choi matrix reveals that the map DTD is completely 
positive if T is. Similarly, if T is merely positive, then DTD is positive as well, since it is a 



concatenation of positive maps. Consequently, the convex combination (T + DTD)/2 = f' inherits 
the property of being (completely) positive from T, ' r-j 

For a qubit map of the form in Eq.([T]i with v = complete positivity can be expressed in 
terms of the singular values of A and its determinant. In order to understand this, suppose that 
one acts with a unitary before and another one after applying the map T so that the overall action 
is p i-» U2T(U\pU\)U2. The matrix representing this concatenation is then given by the prod- 
uct (1 © 2 )f(l © Oi), where O t £ 50(3) are real rotations of the Bloch sphere. This reflects 
the two-to-one group homomorphism SU(2) — > 50(3). Hence, we can use this to diagonalize 
A —> diag(si, s 2 , S3) without changing the complete positivity property. Moreover, we can choose 
the Oi's such that the components of s follow a particular order and satisfy si > s 2 > \ss\. Then 
a generally necessary and for unital maps also sufficient condition for complete positivity is that 
IIVV02I 

si < 1 and si + s 2 < 1 + s 3 . (4) 

These conditions can, alternatively, be summarized by "s £ T" [RSW02] where T C R 3 is the 
tetrahedron whose corners have coordinates with components s, = ±1 satisfying s\s 2 S3 = 1. Note 
that s £ T is independent of the ordering and sign pattern attained from a particular choice of the 
4 G 50(3). 

Since our interest lies, following Eq.pk in the eigenvalues of A but the constraints are essentially 
given in terms of its singular values we need to relate singular values and eigenvalues. This is the 
content of a classic result of Weyl |Wey49| and Horn ||Hor54| of which we need a real version 
provided in [LM01 1. Applied to our context this gives: 

Lemma 1 Given A € C 3 and s £ R 3 , there is a A £ A^R) with spec(A) = A and O1AO2 = 
diag(si, S2, S3) for some Oi € 50(3) iff the following three requirements are fulfilled: 

1. A is, as a multiset, invariant under complex conjugation, i.e., either A G R 3 or there is a single 
real Xi and a complex conjugate pair. 

2 - nLi s * = nLi A «- 

3. \s\s2\ > IAjAjI and \s\\ > |A{| where * refers to a (re-)ordering with decreasingly ordered 



absolute values, i.e., \s^\ > |s^ +1 | and similarly |A^| > |A 



j+ii- 



Proof. Necessity: 1. follows from A being real, 2. follows from both sides being equal to det(A), 
and 3. follows from Weyl's theorem |Wey49| since the |sj| are the singular values of A. 

Sufficiency: Following ILMOTl the conditions 1.-3. guarantee the existence of a matrix M £ 
A^3(R) with spec(M) = A and singular values |sj|. The latter implies that there are Oi £ 50(3) 
s.t. O1MO2 = diag(Si, §2, S3) where |a,-| = |sj|. Moreover, since Yii s i = Yii *i (exploiting 
condition 2.) we have that s and s differ either by no, or by two signs. In both cases we can set 
A := M. If two signs are different we simply multiply one of the 0/s by a diagonal matrix like 
diag(l, —1,-1) or a permutation thereof in order to obtain the O/s. q 

Now we are prepared for solving the inverse eigenvalue problem for qubit channels: 

Theorem 1 (Spectra of qubit channels) Given a multiset A £ C 4 the following statements are 
equivalent: 



1 Note that in higher dimensions the strategy of the proof fails: the higher-dimensional analogue of the above D maps 
p 1— > 2tr [p] t/d — p, which is no longer a positive map if d > 2. Similarly, DTD then can fail to be completely positive 
even if T is so. 



1. There is a trace-preserving completely positive linear map T : A^C) — > A^2(C) such that 
spec(T) = A. 

2. There is a unital and trace-preserving completely positive linear map T : A^2(C) — > A4 2(C) 
such that spec(T) = A. 

3. A = {1}UA where A G C 3 is, as a multiset, closed under complex conjugation. Furthermore, 
if we define s € R 3 by Si :— Xi if Xi G K and Si :— |A.j| otherwise, then (with T the 
tetrahedron defined below Eq.$\) 

s e T. (5) 

Proof. The equivalence 1.-^-2. is a direct consequence of Prop[2] 

3.=*>2.: by Lemma [I] there exists a A G Af 3 (R) with spec(A) = A and OiA0 2 = 
diag(si, S2, S3) for some O, G 50(3). The unital and trace-preserving map defined via T = 1 ® A 
has thus spec(T) = {1} U A = A and is completely positive due to Eq.pl). 

2.=>3.: as discussed above, the matrix representation of the channel has the form T = 1 A and 
spec(T) = {1} U A = A with A = spec(A). Our aim is to argue that diag(l, s l5 s 2 , S3) represents 
a valid quantum channels as well so that Eq.pll has to be fulfilled. Making use of the real Schur 
decomposition we can w.l.o.g. assume that A has one of the following forms 

Ai 

(6) 





where * denotes arbitrary entries and the 'real form' on the left is such that Aj <G R and the 'complex 
form' on the right is such that Ai 6 R and A2.3 = a ± ib with a, b G R. 

Let us consider the real form first. Define a new matrix A' G A^R) via the convex combination 

A ' : = w\ £ DtAD » (7) 

W D t es 

where S is the set of four diagonal matrices D h which have diagonal entries ±1 and satisfy 
det(A) = 1 so that A' = diag(Ai, A 2 , A 3 ). Note that T" := 1 © A' represents a valid quan- 
tum channel since Eq.(|7]l corresponds to a convex combination of channels which merely differ by 
unitary conjugations (represented by 1 Di). Therefore A G T and since in the real case s = A 
Eq.(|5]) follows. 

Consider now the complex form. Again we can get rid of the * -entries by convex combination, if 
we choose S = {l,diag(l, —1,-1)} this time. Since the two-by-two block matrix in the lower right 
corner has singular values | A2 1 , | A3 1 (which equal s 2 and s 3 ) there are special orthogonal matrices 
Oi,C>2 such that OiA'0 2 = diag(si,S2,S3). Consequently, s G T because diag(l, si, S2, S3) 
represents a valid quantum channel since it can be obtained from T by convex combinations and 
concatenating with unitary conjugations. |— I 

A counterpart of this theorem for positive maps can easily be derived: 

Theorem 2 (Spectra of positive qubit maps) . LetT : .A/^C) ~> M.2(C) be a positive and trace- 
preserving linear map. Then the multiset A := spec(T) satisfies 

1. 1 G A and A is invariant under complex conjugation, 

2. \X\ < 1 for all Ag A. 

Conversely, if a multiset A of four complex numbers satisfies these two conditions, then there is a 
unital and trace-preserving positive linear map T on M.i (C) such that spec(T) = A. 



Proof. As discussed before, conditions 1. and 2. are necessary for positive trace-preserving linear 
maps in any dimension. For the converse we construct a map via T := 1 © A with A of the form 
in Eq.d6|, depending on whether or not A = {1, A l7 A2, A3} contains a complex conjugate pair, and 
* = in both cases. By construction we have spec(T) = A and T is trace-preserving and unital. 
Moreover, it is positive since 1 1 A| |oo < 1 guarantees that T maps the Bloch sphere into itself. [— j 



IV. THE NON-ZERO PART OF THE SPECTRUM 

A. The spectral set 

In this subsection we will regard the spectrum as a set. That is, {spec(T)} is a set of complex 
numbers which contains each element only once so that degeneracies are not taken into account. 
This considerably simplifies the inverse eigenvalue problem if we allow for an additional null space: 

Theorem 3 (Spectral sets without zero) Let A = {Ai, . . . , A^} be a set of non-zero complex num- 
bers such that 

1. it is closed under complex conjugation: A = A, 

2. it contains one: 1 € A, and 

3. all elements are contained in the unit disc: |A| < 1, VA G A. 

Then there exists a completely positive and trace-preserving linear map T : Md(C) — > Md(C) 
with d < 2 maxjiV — 1, 1} such that 

{spec(T)}\{0} = A. (8) 

Conversely, if A is the spectral set of a positive and trace-preserving linear map, then conditions 
1.-3. have to be fulfilled. 

Proof. For each Aj G A consider the multiset Aj := {1, 1, Aj, Aj}, According to ThmfT] there is 
a quantum channel Tj : B(rHj) —¥ B(rHj), where Hj ~ C 2 , such that spec(Tj) = Aj. Defining a 
Hilbert space W := (J) -Hj and isometries V 3 ■ : H — > Hj we can construct a quantum channel on 
B(H) via 

T{p) :=Y J V]T i {v jP V})v i . (9) 

3 

By construction Aj e spec(T) for all Aj G A. Moreover, since T (\-ijj) (tfi\} = for all ip E T-Lk, 
cj> G Hi with k j£ I, the kernel completes the spectrum of T such that indeed 

{spec(T)} = AU{0}. (10) 

What is d = dimH? Assume that N r and N c are the numbers of real elements different from 1 and 
complex conjugate pairs in A respectively, i.e., N — N r + 2N C + 1. Then our construction requires 
d = 2(N r + N c ) < 2N - 2 degrees of freedom if N > 2 and d = 2 if A = {1}. 

The converse statement in the theorem is just a restatement of the general properties discussed in 

sec(n] D 



B. Moments and the non-zero part of the spectrum 

From now on we consider the spectrum again as a multiset. That is, every eigenvalue appears 
according to its algebraic multipliticity. A central notion in the following discussion of the non-zero 
part of the spectrum is that of moments. Given a multiset of complex numbers A = {Ai, . . . , A„} 
the fc'th moment of A is defined via 

n 

/i fc (A):=5^Aj, keN. (11) 

j=i 

The collection of moments determines the non-zero part of A and viceversa. Obviously, all moments 

are real if A = A. Conversely, positivity of all moments //£ (A) > 0, k G N implies A = A. Positivity 

also guarantees that A contains a positive element g G A for which VA G A : |A| < g holds l|Fri09l . 

Note that if A is the spectrum of some linear map with matrix representation T, then \x% (A) = 



tr T 



-* 



If T is a linear map on matrices, this can be rewritten as follows: 



Lemma 2 (Moments of maps on matrices) Let T : Md(C) — > A^d(C) be a linear map, then 



tr 



r-ph 



(fi|(T fe ®id)(|fi)(fi|)|0), (12) 



where |fi) := Y] i — 1 \ii). Moreover, ifT is completely positive and has Kraus operators {Kj} then 



tr 



rph 



J2 \tr[K n ---K Jk ] 



(13) 



)!>■■■* 



Proof. Both statements are elementary observations when using matrix units \i) (j | as basis elements 
for computing the trace in Hilbert-Schmidt space as 



,/ 



tr 



rph 



E (*i Tfc (i<)oi)u> 



D 



This leads us to the main theorem of this section: 



Theorem 4 (Non-zero part of the spectrum) Let T : A4d(C) —¥ Md{C) be a trace-preserving 
and completely positive linear map. Then /ife(spec(T)j > for all k G N. 

Conversely, consider a multiset A = {A 1; ...,A„} of non-zero complex numbers such that 
max{|A| A G A} = 1 is attained for a unique element X p G A. If for all k, m G N 

jUfc(A) > 0, and (14) 

Hm(A) > implies /x fem (A) > 0, (15) 

then there exists a primitive trace-preserving completely positive linear map T such that A is the 
non-zero part of the spectrum ofT, i.e., spec(T)\{0} = A. 

Proof. Positivity of the moments follows from Lemma El and the fact that T is assumed to be 
completely positive. For the converse statement we exploit a highly non-trivial result of BBH91I . 
This states that Eqs.( 14|15 l together with the uniqueness of A p are necessary and sufficient for the 
existence of some N £~N and a primitive non-negative matrix M G A4n0&+) withspec(M)\{0} = 
A. From this one can easily construct a stochastic matrix with the same spectrum: denote by {L\ — 



(L\M the left Perron eigenvector and define a diagonal matrix X with entries Xu = (L\i) for a set 
of orthonormal basis vectors {\i)}. Then S := XMX^ 1 is a primitive stochastic matrix which has 
the same spectrum as M. Now define a map T : Mn(C) — > A4jv(C) via 

N 

T(p):= J2 Sij (j\p\j) \i)(i\. (16) 

If S\p) — \\p) , then T (p) — \pfoi p = '^2 li Pi\i){i\ and thus spec(T) D A, i.e., the N eigenvalues of 
S (taking algebraic multiplicities into account) appear in the spectrum of T. Moreover, T ( | i) (j | ) = 
for all i ^ j gives rise to a kernel of dimension N(N — 1) which thus completes the spectrum so 
that indeed A is the entire non-zero part of the spectrum of T. By construction T is trace-preserving 
and completely positive. In fact, it is an entanglement breaking quantum channel. Primitivity of T 
follows from primitivity of S: since S has an entrywise positive fixed point, T has a positive definite 
fixed point. Together with the fact that T has a unique eigenvalue of magnitude one this implies 
primitivity ISP GWC09I . n 

Some remarks on the above theorem are in order: 

1. Positivity does not imply positive moments. Take any positive map T on A^d(C) which is 
not completely positive. Then there is a -0 e C d such that (i^\(T ® id)(|17)(fi|)|?/>) < 0. 
We can always write \ip) — (Y (g) l)|fi) for some Y € A^d(C). Using Eq.(12i the map 

T'{p) := Y^T(p)Y will then be such that tr f" < although T is positive. The same 
holds for n-positive maps with n < d. 

2. The required number of zeros. It was shown in IIBH91I that even for the case n = 4 there 
cannot be an upper bound on the number N — n of zeros which one has to add in order realize 
a given A (which satisfies all necessary requirements) as the non-zero part of the spectrum 
of some N x N non-negative matrix M. A priori, this does not imply that such a bound 
does not exist when the realization is in terms of a quantum channel. In fact, the inequality in 
Eq.( 27 1 from which the necessity of a large number of zeros is usually derived in the classical 
case, fails to hold in the quantum case. And indeed, the additional freedom provided by 
quantum mechanics often allows for much more efficient realizations — see |WPG09| and 
the discussion in Sec. [Vl]for more details. We will, however, see below Thm|6]that also in the 
quantum context ancillary zeros can be necessary in order to realize a certain spectrum. 

In the classical context of non-negative matrices, it has been shown in 1JLL 96 1 that if A is 
the non-zero part of the spectrum of some M <G A4n(WL+) with rank(M) = r, then there is 
anM'e X r 2(R + ) for which spec(M')\{0} = A. That is, the unboundedness of N in the 
classical context is related to large or numerous nilpotent Jordan blocks. 

3. Eq.{15 1 is not necessary. Note that in the converse part of ThmJ4]all the imposed conditions 
on A but Eq.([T5J) are necessary. While for the analogous classical statement for non-negative 
matrices in |BH91 ] Eq.( 15 1 is evidently necessary, too, it fails to be so in the quantum context. 
Let us construct an example of a primitive quantum channel where Eq.([T5| does not hold: 

Define the following d+1 Kraus operators Kj G A4rf(C): 

1 f I j) U + 1 mod d\ for j = 0, . . . d - 1 



I<3 '' v/2 XTfkZlexp (nik/d)\k)(k\ for j = d. (1?) 

The map T(p) := J2j=qKjpK'- is then completely positive, trace-preserving and unital. 



Moreover, following Eq.( 13 i we get tr 



T 



if d > 2. Primitivity of 



> although tr 
the map T is guaranteed since for n > rfwe have 

spanji^ ■■■K jn } = M d (C) (18) 



which implies primitivity by virtue of [SPGWC09|. 



Although tr 



T ; 



in the above example, we note that primitivity generally implies that 



tr 



rpk 



> for all k > p where p depends on the dimension only. This follows from 

Wielandt's inequality which in the classical case of stochastic matrices in MdO&+) yields 
p — d 2 — 2d+2 and in the context of quantum channels on A^C) it gives p = (d 2 — 2)(d 2 — 1) 
HSPGWC09L 

4. Classical vs. quantum channels. Consider a matrix representation T of a map T : A4 d(C) — >■ 
A4d(C) expressed in terms of matrix units, i.e., (k,l\T\i,j) := (l\T(\i)(j\)\k) for i,j, k,l = 
l,...,d. Then the d x d matrix S defined via Sij := (ii\T\jj) is entry-wise non-negative 
if T is a positive map and it is stochastic if in addition T is trace-preserving. That is, every 



quantum channel T contains a stochastic matrix S as a principal submatrix. Eq.( 16 1 reverses 
this: it defines a quantum channel from this stochastic submatrix by setting all other matrix 
elements of T to zero. 

The close relation between the non-zero spectra of quantum channels and of stochastic matrices 



implied by Eq.(16i allows to exploit other results from the classical world. In the following we 
derive in this way two sufficient conditions for realizable spectra for which a bound on the number 
of ancillary zeros can be given: 

Theorem 5 (Non-zero spectrum with non-positive real parts) Let A = {Ai, . . . , A n } be a multi- 
set of non-zero complex numbers with Ai = 1 and Rc(Aj) < 0/or all j > 2, and let d € N be the 
smallest number such that /iti(A) 2 < d/i2(A). If A = A and 

Mi(A)>0 and p, 2 {A) > 0, (19) 

then there exists a primitive, completely positive and trace-preserving linear map T : A^^C) — > 
Md(C) such that A = spec(T)\{0}. 



Proof. The theorem is a direct consequence of its classical counterpart derived in |LS06| inserted 
into Eq.(jT6]>. As discussed already in the proof of ThmHIprimitivity is carried over from the classical 

case. [-] 

Theorem 6 (Real non-zero spectrum for n = 4) Let A = {Ai, . . . , A4} be a multiset of non-zero 
real numbers. If (i) A = A, (ii) for all j : \Xj\ < 1 € A and (Hi) X)i=i A; — 0> tnen tnere 
exists a completely positive and trace-preserving linear map T : A^C) — > A'f^C) such that 
A = spec(T)\{0}. 

Proof. Following IILL781 the stated conditions are sufficient for the existence of a non-negative (and 



by similarity transformation stochastic) matrix with spectrum A. Using Eq.( 16 1 again this gives the 
desired result. |— I 

Thm|6]together with ThmfTlshows that ancillary zeros can be necessary to make a certain spectrum 
realizable. Consider for instance A = {1,1,1,-1}. Then ThmfTI implies that this cannot be the 
spectrum of a quantum channel on A^C) whereas Thm|6] tells us that it can be the non-zero part 
of the spectrum of a quantum channel on A^C). 

Finally, we use the similarity between the classical and the quantum case to characterize the non- 
zero part of the spectrum of quantum channels with 'full Kraus rank'. The latter means that the 
complex linear span of the Kraus operators of a quantum channel on M. d {£.) is the entire matrix 
space A'fd(C). Equivalently, the Choi matrix (or Jamiolkowski state) has full rank. This condition 
can be seen as the quantum counterpart to entry-wise positivity of stochastic matrices. Note that, in 
particular, every quantum channel with full Kraus rank is primitive. 
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Theorem 7 (Non-zero spectrum under full Kraus rank) Let A = {1,\ 2 . . . X n } be a multiset of 
non-zero complex numbers. There is a completely positive and trace-preserving linear map T with 
full Kraus rank such that A = spec(T)\{0} iff \ Xj \ < I for all j ' > 2 and pk (A) > Ofor all k G N. 



Proof. Necessity of the conditions is implied by Eq.( 12 1 and the general discussion in Sec|ll] 



Following |BH91| the stated conditions are sufficient for the existence of an entry- wise positive 
matrix M with A = spec(M)\{0}. As mentioned in the proof of Thmplthis can, by similarity 
transformation, be made a stochastic matrix S which is again entry-wise positive. Inserting this 
into Eq.([T6} the latter guarantees full Kraus rank since the Choi matrix of the constructed quantum 
channel reads 

(T®id)(|nKn|)=^s y |i)(i|®|j)<j|>Bi, (20) 

ij 

where s = min{ Si j} > 0. rj 

V. CYCLES AND THE PERIPHERAL SPECTRUM 

We will now have a closer look at the space corresponding to the peripheral spectrum (i.e., eigen- 
values which are phases) of a positive linear map T : Md(C) — > At<j(C) which is assumed to be 
either trace-preserving or unital. Denote the complex linear span of all the respective eigenspaces 
by 

X T := spanJAT G M d (C)\3ip e R : T(X) = e iv X), (21) 

for which we will simply write X if the dependence on T is clear from the context. The questions 
which we address first are: what's the structure of XI what's the structure of T restricted to XI and 
what does this imply for the peripheral spectrum? 

Recall that if T is trace-preserving and positive, there cannot be any non-trivial Jordan-block 
associated to an eigenvalue of modulus one. 

The following shows that X is a *-algebra (w.r.t. a modified product) and that T acts on it either 
by permuting blocks or by unitary conjugation on subsystems. The property which we use in order 
to prove the assertion lies in between positivity and complete positivity: we call a positive trace- 
preserving map T on M d (C) a Schwarz map if T*(A^A) > T*(A^)T*(A) holds for all A € 
■Mrf(C). This is in particular true if T is completely positive. 2 

Theorem 8 (Structure of cycles) Let T : A^d(C) — » A^d(C) be a trace-preserving Schwarz map. 

1. There exists a decomposition of the Hilbert space C d — Ho © ®fc=i ^fc ^ nto a direct sum of 
tensor products^ H). = Hk,\ ® T~Lk.2 and positive definite density matrices pk acting on 1ik,2 
such that 

K 
X T = 0®($M dk ® Pk , (22) 

fe=i 



2 In fact, 2-positivity, i.e., positivity of T Cg) id2, is sufficient but not necessary for T to be a Schwarz map. 

3 Note that we may have dirnHo = as well as dimHk.i = 1. Moreover, the direct sum in the decomposition does not 
necessarily correspond to a block structure in computational basis -just in some basis. 
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where Md k is a full complex matrix algebra on %k,i of dimension dk '■= dim('Hfc.i)- That is, 
for every X G Xt there are Xk G B{'Hk,i) such that 

I = O®0x fc ®pi. (23) 



2. There exist unitaries Uk G B(Hk,i) and a permutation n, which permutes within subsets 
of {1, . . . , K} for which the corresponding Hk.i 's have equal dimension, so that for every 
X G Xt represented as in Eq. ( |23| ) 

K 

r(i) = oe0(/ Aw (/|® ft . (24) 

fc=i 

Proof The basic ingredients of the proof are those which already appeared in [WC08 1 where pos- 
itive maps with unit determinant were characterized. Exploiting Dirichlet's lemma on Diophantine 
approximations 4 we can find an ascending subsequence rij G N so that lim.j_). 00 T ni converges to 
a map / with eigenvalues which are either one or zero. Clearly, / is again a trace-preserving, posi- 
tive linear map which satisfies the Schwarz inequality since all these properties are preserved under 
concatenation. Moreover, Xt = Xj is the fixed-point space of I. So, statement 1 . of the theorem 
follows from the structure of the space of fixed points provided in MLin99ll . 

2. The map T _1 := rim,_ Kx) T" i_1 is the inverse of T on X since by construction T~ l T = I. 
Hence, both T and T _1 map X — >• X in a bijective way. The crucial point here is that T _1 is a again 
a positive, trace-preserving map since it is constructed as a limit of such maps. To understand the 
consequences, consider any pure state in X, i.e., a density operator a G X which has no non-trivial 
convex decomposition within X. Then the image of a under T has to be a pure state as well: assume 
this is not the case, i.e., T(a) — £\ ^i a i is a non-trivial convex decomposition into states at G X. 
Then applying T _1 to this equation leads to a contradiction since a — J2i AjT 1- (<7j) is not pure. 
Consequently, both T and T _1 map pure states in X onto pure states. Note that a pure state can 
only have support in one of the K blocks, for instance a — x p^ where x G B(Hk,i) is a rank 
one projection. Now we know that T(a) = x' ® pt' for some k' and some rank one projection 
x' G B(Hk'.i)- By continuity and the fact that T is a bijective linear map on X, we have that within 
X: (i) every element of block k is mapped to an element of the block k' (i.e., k' does not depend on 
x), and (ii) the spaces Hk,i and Hk/.i must have equal dimension. Therefore, there is a permutation 
7r which permutes blocks with equal dk so that X G X is mapped to 

r(I) = O90T ft (a; l(l) )ep S) 

k 

with some linear maps Tk : M.d k ~ > Md k - Since the latter have, together with their inverses, 
to be positive and trace-preserving they must be either matrix transpositions or unitary conjuga- 
tions |WC08 1. Matrix transpositions are, however, ruled out by the requirement that T and thus each 
Tk is a Schwarz-map. |— I 



The action of T on Xt derived in Eq.(24 1 now determines the possible peripheral spectra: 



Theorem 9 (Peripheral spectrum) Let T : Md(C) —> Md{C) be a trace-preserving Schwarz 
map. There are integers n c ,d c G N satisfying ^2 c n c d c < d, integers m c G Z„ c and vectors 



4 Dirichlet's lemma states the following |Sch80|: let x £ R m and q > 1 any integer. Then there exist integers 
n, pi , . . . , pm such that 1 < n < q m and \xkn — Pk\ < !/</ for all k. 
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y«' c ' G C dc w/iose components are phases (i.e., Vfc : |/i). | = 1) so ?/m£ the peripheral spectrum of 
T (including algebraic multiplicity) is given by the ^ c n c d 2 c numbers 

(c) -(c) 2 " m <= , OC s 

Mfc Mj e "= . (25) 

Conversely, given such a set of numbers there is a completely positive, trace-preserving and unital 
linear map T on M.d(C) with d = J2 C n c d c so that its peripheral spectrum coincides with the given 
numbers and it has no other non-zero eigenvalue. 

Proof. First note that the eigenvalue equation T(X) = XX for X G Xt is, via Thm|8] equivalent 
to 

Vfc : x <k) = XUlx k U k . (26) 

Assume for the moment that 7r is a cycle of length n c , i.e., after relabeling we have fc <E Z„ c and 
7r(fc) = fc + 1 mod n c - the general case will later reduce to this one by decomposing a general 
permutation into cycles. All C/^'s and x^'s now act on spaces of equal dimension, d c say. For given 
ZTfc's and A the eigenvalue equation determines each x k from xq and leads to the additional constraint 
x Q = \ n "U^x U with U :=U U 1 ---U ne -i. Now define any unitary U on C dc for which U n " = U 
and consider any of the d 2 . solutions of the eigenvalue equation y = jiWyU, y 6 B(C dc ). Then 



xq = y and A = /i cxp [2i:im/n c ] leads to a solution of Eq.( 26 1 for every m 6 Z„ c . By construction 
the eigenvectors X of T which correspond to different m's will be orthogonal, so that in total the 
construction leads n c d 2 c and thus all solutions. If /ii, . . . , /id c are the eigenvalues of U, then the 
multiset of eigenvalues A is given by {^fcMi ex P [27rim/n c ]} with fc, j = 1, . .. ,d c and m G Z„ c . 

The case of a general permutation now follows by decomposing it into cycles (labeled by c and 
of possibly different lengths n c , which implies that ^ c n c — K) and observing that the peripheral 
spectrum of T is just the union of the spectra of all the 'cycle-building-blocks' . 

In order to prove the converse we construct a map T guided by the structure appearing in Thm|8] 
(but now with dim - Hfc2 = 1). The desired T is obtained as a concatenation of three types of maps: 
a pinching onto the ^ c n c d c blocks of C d = (J) c C rfc ®" c , a permutation with cycles of length n c 
which permute blocks of dimension d c and, finally, unitary conjugations on each block. In this way 
we achieve that T is (i) completely positive, unital and trace-preserving, (ii) it projects .M d(C) onto 
a direct sum of matrix algebras of the form in Eq.([22]i (albeit with pk = 1 and without kernel) 



and (iii) it acts on this sub-algebra as in Eq.(24 1. It follows then from the previous analysis that an 
appropriate choice of the unitary conjugations will lead to the predetermined spectrum. [— ] 

Note that a different converse can be proven along the same lines as in Thm|9] for every trace- 
preserving Schwarz-map T there is a trace-preserving completely positive map T such that T = T 
on Xt- 5 That is, complete positivity is not only irrelevant for the peripheral spectrum but for the 
entire peripheral eigen-system. As shown below ThmHlthis does no longer hold for the rest of the 
spectrum, let alone the corresponding eigen-system. 



VI. CLASSICAL VS. QUANTUM INVERSE EIGENVALUE PROBLEMS 

We will briefly review the classical counterpart of our treaties — the inverse eigenvalue problem 
for stochastic matrices to which we resorted already in Sec |IVB| For a more detailed overview see 
BCG05L Recall that the spectrum of any stochastic matrix is contained in the unit disc, it contains 



5 This is closely related to the well known fact that any conditional expectation, i.e., positive projection onto a *-subalgebra 
of A1d(C), is automatically completely positive. 
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one, and it is, as a multiset, invariant under complex conjugation. A non-negative matrix 6 is ir- 
reducible iff the spectral radius is a non-degenerate eigenvalue and the corresponding eigenvector 
is entrywise positive definite. An irreducible matrix in turn is primitive iff there is only a single 
eigenvalue whose modulus matches the spectral radius. 



A. The classical inverse eigenvalue problem 

Every stochastic matrix S is, after a permutation, upper block-triangular with the diagonal blocks 
being irreducible. That is, the spectrum of S is the union of spectra of irreducible substochastic 
matrices. As used in the proof of ThmH] already, every irreducible matrix is, up to a rescaling, 
similar to a stochastic matrix. For every irreducible matrix M € A4 ( j(R+) there is an index of 
cyclicity k € {1, . . . , d} so that the spectrum of M is invariant under 27r/fc-rotations in the complex 
plane. That is, the eigenvalues come in k rotated d/fc-tuples. Correspondingly, M k is block-diagonal 
with primitive matrices in the diagonal blocks. So finally, the spectrum of any non-negative matrix, 
stochastic or not, is essentially determined by the spectra of primitive matrices. 

The non-zero part of the spectrum of primitive matrices is completely characterized by Eqs.( 14|15 1 



together with the fact that the spectral radius is attained only for a single eigenvalue BBH91L The 
number of additionally required zero eigenvalues is, however, generally unbounded. One way to 
show this is to exploit the dimension dependent condition 

Mfe(A) m < d m - Vfcm(A), k, m e N (27) 

which is valid for any non-negative matrix in A4d0&+) IILL781lJoh81l . For a multiset A of three 
complex numbers Eqs.([T4J27i are in fact necessary and sufficient for A being the spectrum of a 



matrix in A^3(R+) |LL78|. For d — 4, however, these conditions are no longer sufficient. 



B. Comparing classical and quantum 

Maybe the most basic difference between the spectra of classical and quantum channels concerns 
the location of single eigenvalues |WPG09|. While for quantum channels on A / (d(C) the entire 
unit disc is accessible for any d > 2, stochastic matrices in A4d(Si+) obey a dimension dependent 
constraint on the location of every single eigenvalue: a necessary condition for A to be a valid 
eigenvalue is that it is contained in the convex hull of all roots of unity of order up to d. A complete 
characterization of the accessible region within the unit disc can be found in [Kar51 DD46]. 

As a consequence, for certain sequences of four non-zero eigenvalues which obey Eq.([T4| there 
is a quantum channel on Ai 2 ('C) but no classical channel for any bounded dimension. That is, while 
on the classical side an unbounded number of zeros have to be added, no extra kernel is required on 
the quantum side. 

As discussed below Thm|4] already, it is, however, open whether or not there is a general bound 



on the number of required zeros in the quantum case. Eq.(27i, which is classically used to bring 



the dimension into play, fails to hold for spectra of quantum channels: the channel constructed in 



Eq.( 17 1 leads to /ii > with /j,2 = and thus violates Eq.( 27 1 for k = 1, m = 2. 

Another crucial difference between the classical and quantum case is the breakdown of the prob- 
lem into irreducible and finally into primitive blocks. The decomposition valid for stochastic ma- 
trices fails to hold for quantum channels and no counterpart of this is known for the quantum case. 
Maybe the structure of the peripheral eigen-system proven in Thm|8]hints at a more general decom- 
position. 



1 A matrix is said to be 'non-negative' if all its entries are non-negative. 
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VII. APPLICATIONS FOR TIME SERIES AND CORRELATION FUNCTIONS 

Suppose we are given a bounded sequence a £ l°° which is generated via 

at = {A\f t \p), teM , (28) 

from powers of a linear operator T e M.r>{C). This may for instance be the discrete time evolution 
of some expectation value or, as explained later, a two-point correlation function. Then a contains 
information about the spectrum of T. A simple relation between spec(T) and a is provided by the 
Cayley-Hamilton theorem, i.e., the fact that T is a solution of its characteristic equation. So let 

p( x ) = IIfc=i(^fc~ x ) — Sfc=o Ckxk b e the characteristic polynomial ofT, then J2k=o a k+i c k = 
for all I 6 No. Depending on degeneracies of the eigenvalues and their Jordan-block structure there 
can be annihilating polynomials of degree smaller than D IGa n74al . 

Conversely, eigenvalues of T can be determined as roots of a polynomial which is constructed 
directly from the sequence a. To this end define for each r <G No a vector v T 6 l°° via v T := 
(a T , a r +i, . . .). If the real linear space V := span{?; T }i-eN has dimension r < oo, then there are 
coefficients Ck such that for all I £ No'- vi +r — Y^k=o CkVl + k - These coefficients, which can be 
determined from a, define a polynomial p(x) := x r — X)I=o Ckxk '• The roots of p then turn out 
to be eigenvalues of T (since they coincide with the poles of the z-transform discussed below, see 
HGan74bl p. 205 ff). 

A related connection between spec(T) and a is provided by the z-transform C : C — > C: 

c ^--=\T. a i - (4(*i-*r>>. (29) 

tGN 

Note that the series which defines C converges outside a disc with radius equal to the spectral radius 
of T and C is defined inside by analytic continuation. 

The poles of C correspond to eigenvalues of T. Depending on the interplay between T, A and p, 
not all eigenvalues may, however, appear as poles. Consider as an example the sequence 



a t = (lm[(l-i) t x t ]+y t )/V2, x, y£ [-1,1]. (30) 



In this case C has poles at zi t % = (1 ± i)x and at z?, = y. So according to ThmfTJthe sequence can 
be generated by a quantum channel acting on A / J2(C) only if (\z\\, \z2\, 23) € I. In fact, it can be 
generated by a map on M.%(C) with matrix representation 

f=( 1 °W x x 

V ° y J \-x x 

with p and A given by (1, 1/v2j l/\/2,0) and (0, 1,0, 1) respectively. So the example can be 
generated by a qubit channel iff the poles are valid eigenvalues. Also note that T is a positive map 
iff \y\ < land |*| < 1/V2. 

Using the above relations, knowledge about the spectra of quantum channels can be used to decide 



whether or not a given sequence a has a representation of the form in Eq.( 28 1 with T being a quantum 
channel of a certain dimension. Such sequences arise naturally for instance as discrete, homogeneous 
time-evolutions. They also appear as two-point correlation functions in quantum spin chains whose 
state is finitely correlated |FNW92 1. The Kraus operators of the T then correspond to the matrices of 
the matrix product state |PGVWC07| and t + 1 is the distance between the two considered sites on 
the chain. So a detailed spectral analysis might help in this context if one wants to find an accurate 
finitely correlated approximation to given two-point correlations. 
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FIG. 1: Sequences given by Eq.l|30| with y — 2/3. While the solid curve (x — 0.58) can be generated by an 
evolution of a qubit quantum channel, the dashed curve (x — 0.59) cannot since the poles of its z-transform do 
not correspond to valid eigenvalues of a qubit quantum channel. The sequence can, however, be generated by a 
positive trace-preserving map. 
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